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ABSTRACT 

The relationship between galaxies and dark matter can be characterized by the halo 
mass of the central galaxy and the fraction of galaxies that are satellites. Here we 
present observational constraints from the SDSS on these quantities as a function of 
r-band luminosity and stellar mass using galaxy-galaxy weak lensing, with a total 
of 351 507 lenses. We use stellar masses derived from spectroscopy and virial halo 
masses derived from weak gravitational lensing to determine the efficiency with which 
baryons in the halo of the central galaxy have been converted into stars. We find that 
an galaxy with a stellar mass of 6 x 10^°Mq is hosted by a halo with mass of 
1.4 X lO^^/i~^M0, independent of morphology, yielding baryon conversion efficiencies 
of ITlg^ (early types) and 16^g^ (late types) per cent at the 95 per cent CL (statisti- 
cal, not including systematic uncertainty due to assumption of a universal initial mass 
function, or IMF). We find that for a given stellar mass, the halo mass is independent 
of morphology below Msteiiar = IO^^Mq, in contrast to typically a factor of two dif- 
ference in halo mass between ellipticals and spirals at a fixed luminosity. This suggests 
that stellar mass is a good proxy for halo mass in this range and should be used pref- 
erentially whenever a halo mass selected sample is needed. For higher stellar masses, 
the conversion efficiency is a declining function of stellar mass, and the differences in 
halo mass between early and late types become larger, reflecting the fact that most 
group and cluster halos with masses above IO-'^^Mq host ellipticals at the center, while 
even the brightest central spirals are hosted by halos of mass below IO^^Mq. We find 
that the fraction of spirals that are satellites is roughly 10-15 per cent independent 
of stellar mass or luminosity, while for ellipticals this fraction decreases with stellar 
mass from 50 per cent at IO^^Mq to 10 per cent at 3 x IO^^Mq or 20 per cent at the 
maximum luminosity considered. We split the elliptical sample by local density, and 
find that at a given luminosity there is no difference in the signal on scales below 100 
/i~^kpc between high and low density regions, suggesting that tidal stripping inside 
large halos does not remove most of the dark matter from the early type satellites. 
This result is dominated by halos in the mass range 10^'^ — lO^^ft-^^M©, and is an 
average over all separations from the group or cluster center. 

Key words: galaxies: haloes - galaxies: stellar content - gravitational lensing. 



1 INTRODUCTION 

The connection between the spatial distribution of galax- 
ies and dark matter (DM) is an essential ingredient in 
* Electronic address: rmandelbOprinceton.edu the physics of galaxy formation. One very useful probe of 
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the galaxy-DM connection that recently became available 
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tational lensing induces tangential shear distortions of back- 
ground galaxies around foreground galaxies, allowing di- 
rect measurement of the galaxy-DM correlation function 
around galaxies. The individual distortions are small (of 
order 0.1%), but by averaging over all foreground galax- 
ies within a given subsample, we obtain high signal to noise 
in the shear as a function of angular separation from the 
galaxy. If we know the lens redshifts, the shear signal can be 
related to the projected mass density as a function of proper 
distance from the galaxy. This allows us to determine the 
averaged DM distribution around any given galaxy sample. 

In recent years, the progress on the observational side 
of g-g lensing has been remarkable. In the latest Sloan 
Digital Sky Surv ey (SDSS) analyses JSheldon et alJ 12004 
ISeliak et al Eooi)," 20-30(7 detections of the signal as a func- 
tion of physical separation have been obtained. Similarly 
high S/N detections have also been observe d as a function 
of an gular separation with other surveys jHoekstra et alJ 
|2004) , but the ability to use spectroscopic redshifts for lenses 
is a major advantage to doing lensing with the SDSS. The 
high statistical power has been accompanied by a more care- 
ful investigation of systematic errors, such as calibration bi- 
ases and intrinsic alignments, which for the SDSS are cur- 
rently around 10 per cent and therefore al ready dominate 
the error budget iMandelbaum et alJl2005ah . 

In this work, we seek to use g-g weak lensing to ex- 
plore the galaxy-DM connection for particular subsamples 
of lenses. By comparison with the predicted signa l from a 
halo m odel, as done for simulations in iMa ndclba um et alJ 
JioolH), we can extract average central halo masses and 
satellite fractions. These calculations are done as a func- 
tion of morphology and of environment, in samples se- 
lected based on both stellar masses and luminosity. We 
expect that the divisions by morphology and by envi- 
ronment may be related, due to the relationships be- 
tween color and environme nt, with red galaxies typically 
found in overdense regio ns jDavis fc Gelleil 119761: iDressleJ 
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By determining average central halo masses and satellite 
fractions as a function of these parameters, we hope to gain 
some insight into processes of galaxy formation, and ulti- 
mately into the galaxy-DM connect ion. We note that this 
approach has bee n used before , bv [ Guzik fc SeliakI ll2002ll 
based on data in iMcKav et al ] J200l[) . but with a simpler 
form of the halo model, with a much smaller sample of lenses 
so lower statistical power, and only using the luminosities, 
not stellar masses. Due to our larger sample of lenses, our 
better-understood calibration, and our inclusions of stellar 
masses which are better tracers of stellar and dark mat- 
ter content than luminosities, this work constitutes a sig- 
nificant improv e ment over that one. Another recent work, 
iHoekstra et alJ l|200^, used stellar masses for RCS data 
derived from B — V colors from CFHT photometry in or- 
der to derive halo masses as a function of luminosity, and 



star formation efhciencies as a function of morphology, but 
used only isolated lenses and thus did not derive satellite 
fractions. Furthermore, the lack of spectroscopic redshifts 
for lenses in that work, which allow the derivation of stel- 
lar masses via spectral indicators and the computation of 
the lensing signal as a function of transverse separation in 
this work, complicates the analysis. Halo model analysis of 
galaxy-galaxy autocorrelations has been done obsorvation- 
ally by several groups (.van den Bosch et al. 2003: Cooravl 
I2OO5I IZehavi et al J 12005^ 7 and this halo model analysis of 
galaxy-DM cross-correlations is, in many ways, complemen- 
tary to that approach. 

We begin by introducing the g-g lensing formalism and 
the halo model that is used to extract information about 
central halo masses and satellite fractions in ^ 331 includes 
a description of the SDSS data used for this analysis. We 
present the lensing signal and the halo model fits in ^ and 
interpretation of these results. We conclude in JS] with a 
summary of our findings. 

Here we note the cosmological model and units used in 
this paper. All computations assume a flat ACDM universe 
with Qrn ~ 0.3, f^A = 0.7, and as = 0.9. Distances quoted 
for transverse lens-source separation are comoving (rather 
than physical) /i~^kpc, where Ho = 100 /i kms~^Mpc~^. 
Likewise, AS is computed using the expression for Tj^^ in 
comoving coordinates, Eq.0] In the units used. Ho scales out 
of everything, so our results are independent of this quantity. 
All confidence intervals in the text and tables are 95 per cent 
confidence level {2a) unless explicitly noted otherwise. 



2 WEAK LENSING FORMALISM AND HALO 
MODEL 

Galaxy-galaxy weak lensing provides a simple way to probe 
the connection between galaxies and matter via their cross- 
correlation function 



(1) 



where Sg and 5m are overdensities of galaxies and matter, 
respectively. This cross-correlation can be related to the pro- 
jected surface density 

E(J?) =p J [l + ^g,m [VR^ + X')] dx (2) 

(where = R^ + x^), which is then related to the observable 
quantity for lensing, 

AE(i?)=7t(7?)E, = E(<7?)-S(7?), (3) 

where the second relation is true only for a matter distri- 
bution that is axisymmetric along the line of sight. This 
observable quantity can be expressed as the product of two 
factors, a tangential shear -yt and a geometric factor 

_ Ds , . 

" 4-kG DLDLsil + ZL^ ^ ' 

where Dl and Ds are angular diameter distances to the lens 
and source, D^s is the angular diameter distance between 
the lens and source, and the factor of (1 + zl)~'^ arises due 
to our use of comoving coordinates. For a given lens redshift, 
rises from zero at Zs = z_l to an asymptotic value at 
Z3 3> zl; that asymptotic value is an increasing function of 
lens redshift. 
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There are two basic approaches that can be used to ex- 
tract information about properties of the galaxy distribution 
(e.g., halo masses and satellite fractions) from the g-g lens- 
ing signal. The first approach is to compare directly against 
N-body simulations JGuzik fc Seliakl200ll:lYang et aljl2003t 
iTasitsiomi et alJl2004l:IWeinberg et al.ll2004l) . While this ap- 
proach has the advantage of being fairly direct, it has the dis- 
advantage that even assuming a given cosmological model, 
the process of galaxy formation is not sufficiently understood 
to result in unique predictions for that model. When com- 
bined with the fact that the cosmological model is not itself 
fully determined yet, one would have to expend tremendous 
computational resources to run multiple simulations with 
different cosmologies and models. Furthermore, current sim- 
ulations still suffer from a limited dynamical range, in the 
sense that they require a high mass and force resolution to 
resolve individual galaxies and their associated DM halos, 
while at the same time they must also have sufficiently large 
volume to simulate a representative region of the universe. 
Several simulations of varying box size are thus needed to 
cover the whole observational range in luminosity and scale. 

Another approach to model the relation between 
galaxies and dar k mat t er is to use a halo model (e.g. 
Peacock fc SmithI 120001: ISeliakI I2OOOI : IScoccimarro et alJ 



2OOII: ICoorav fc Shethll2002^ . which in this application can 



be used as a phenomenological description of the processes 
that determine the lensing signal for particular types of 
galaxies. Comparison of the halo model predicted signal ver- 
sus the real signal can lead to the determination of quan- 
tities such as the virial mass distribution and the fraction 
of these galaxies that are satellites, which are useful quan- 
tities for constraining the galaxy formation models and cos- 
mological models. This approach has the advantage that it 
is not as computationally expensive, so large areas of pa- 
rameter space can be explored very quickly, but must be 
com pared against simulations to ensure that it works prop- 
erly. iMandelbaum et alj ll2005lJ) compared the halo model 
with simulations, and determined that central halo masses 
(A/cent) and satellite fractions {a) can be extracted ade- 
quately from the lensing signal, provided that the distribu- 
tion of central halo masses is not too broad (ctm jM^mt less 
than about a factor of two), or that corrections be applied 
if it is broad. Thus, in this paper we use those results to 
extract information from the measured weak lensing signal 
from th e SPSS using halo model fits, as has already been 
done by ISeliak et alJ |2005) (including variation of the re- 
sults with 0.m and erg, for which we refer the reader to that 
paper). The halo model us ed here is the same as that from 
iMandclbaum et al.l (l2005d) . to which the interested reader 
can refer for details. Here we give only the basic details nec- 
essary to understand this paper. 

The halo model can be used to derive the galaxy-DM 
cross-power spectrum Pg,dm{k) by considering separately 
the contributions of central galaxies (which lie in halos 
that are not contained within another halo) and satellite 
galaxies (which lie in subhalos contained entirely within 
another halo). There are then one-halo or Poisson terms, 
derived from the cross-correlation of the galaxy with its 
own matter distribution (for central galaxies and satellites) 
and with that of the host halo (for satellites), and a halo- 
halo (h-h) term derived from the cross-correlation with the 
mass of other halos. The latter is negligible on the small 



(< 2 ^~^Mpc) scales to which we limit ourselves in this work. 
Therefore, rather than attempting to fit for it by fitting to 
the bias, we fix the bias by doing the halo model fits with- 
out the h-h term t o get the halo mass, then using b{M) 
from lSeliak fc Warr en (2004) to redo the fits for Mcent and 
a with the appropriate value of bias fixing the h-h term. 
Including this correction makes very little difference on the 
final results, as expected since on small scales the h-h term 
is negligible. For the same reason we ignore the morphology 
dependence of h alo bias, suggested by recent simulations 
iGao et all ll2005l) that studied the dependence of clustering 
on halo formation time. For splits by local density, we use the 
appropriate b{Ad) for that mass for the high-density samples 
(again, ignoring dependence of bias on local environment), 
and do not include a h-h term at all for low-density samples. 

The one-halo term requires various ingredients, such as 
the halo mass function dn/dM, the radial profile of dark 
matter within halos p(r), the radial distribution of galax- 
ies within groups and clusters, the conditional halo mass 
probability distribution, and the tidal stripping of satellites 
in clusters. T he DM pr ofiles are assumed here to be NFW 
profiles (Nava ,rro^^|T 996'l with co ncentration parameter 
feuUock et al.ll200ll:lEk"e'et alJl200lll 
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where Mni(z), the nonlinear mass scale, is defined such 
that the rms linear density fluctuation extrapolated to red- 
shift z within a sphere containing mass M^i is equal to 
5c — 1.686, the linear overdensity at which a spherical per- 
turbation collapses. Mni(«) is a cosmology-dependent pa- 
rameter. The satellites are assumed to be distributed ac- 
cording to an NFW profile with concentration parameter 
Cg = Cdm ('Carlberg et al' ^199% iPeacock fc S mith' '2000|; 
iGao et al. 2004a) , though the im portance of this a ssumption 
- which may overestimate Cp ||CarIber£,et al. 199j[: lLin et alJ 
120041: iHansen et allbOOSi: iNaeai fc Kravtsov..2005ll) - wiU be 
explored later. Unfortunately, since Cg is, in general, poorly 
determined from previous data and our own, we cannot fit 
for it and must assume some model. 

The conditional halo mass probability distribution for 
a given lens luminosity Li is modeled as having two parts, 
p'^ (central) and p^*^ (non-central), using a free parameter 
Q, the satellite fraction: 



piM■L^) = (1 



a)p'^{M-Lr)+ap^'^{M-L,). (6) 



For halos hosting a central galaxy, we model the luminosity- 
halo mass relationship as a delta-function; for halos host- 
ing non-central galaxies, we assume that there is a relation- 
ship between the number of galaxies of this luminosity and 
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we use {N{M)) cx M' with e = 1 for Af > 3Mcent and 
e — 2 below that value. T his result was sh own there to 
match the simulations from ITasitsiomi et all (j£o04) quite 
well, and the power-law exponent e = 1 above a cut-off 
is consistent with semi-analytic models of galaxy formation 
(Kauffmann c t alJll99gl:lGuzik fc S ehS" 20021. N-body sim- 
ulations (iKravtsov et al.l |200^JZehjvi et^aL, 2005^'l. and ob- 
servational measurements Jjine et alJll998l : IlIu et alj|2004l: 
ICooravl2005l:IZehavi et alJl2005ll . Since some of the observa- 
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tional analyses (e.g. lZehavi et ailbOOST) suggest that signifi- 
cantly different values of e may still be allowed, we explore 
in this paper the dependence of our results on its assumed 
value. 

Finally, our model for tidal stripping is that satellites 
have half of their mass stripped, which corresponds to trun- 
cating the central density profile at 0.4r„ir (beyond which 
we use AE oc as for a point mass). This assump- 

tion is consistent with the average mas s loss for subha- 
los observed in cosmolq gical simulations llGao et alJl2004bt 
iNaeai fc Kravtsovll2005l) . We will attempt to say something 
about this assumption by looking at galaxies in low-density 
and high-density regions separately. 

After using these inputs to obtain the galaxy-DM cross- 
power spectrum Pg.dm{k), we can Fourier transform to ob- 
tain the correlation function S,g^dm{R), integrate once to get 
E(7?) via Eq. |5| and integrate again to obtain AE(_R) via 
Eq.Ol 

Using the precomputed signal, we compare against the 
lensing signal to derive two properties of the lens sample: 
the halo mass Alcent of the central galaxy, defined in terms 
of the radius such that the overdensity within is equal to 180 
times the mean density (roughly 30 per cent larger than the 
mass M200 defined within the radius where the overdensity 
is equal to 200 times critical density, another definition of 
virial radius and mass that is often used); and the fraction 
of galaxies in the sample that are satellites residing in a 
larger halo (i.e., a group or cluster), a. Note that we do not 
attempt to measure any other parameters of the central halo 
mass distribution besides the central halo mass. In these fits, 
the average redshift of the sample is relatively unimportant 
compared to the cosmology {Qm and erg), for wh ich variation 
of be st-fit satellite fraction may be significant iSeliak et alJ 

liooi). 

One issue raised in iMandelbaum et alJ i2005lJ) is the 

question of the meaning of the best-fit halo mass; for sam- 
ples of lenses with broad central halo mass distributions, 
such that the mean and median are significantly different, 
the best-fit masses must be adjusted upwards if one is in- 
terested in the mean (as we are), and downwards if one 
is interested in the median. We attempt to isolate samples 
of lenses with narrow distributions in central halo mass in 
two ways: first, by using bins narrow in luminosity (half or 
one absolute magnitude or a factor of 1.6-2.5 in luminosity) 
or stellar mass (a factor of 2 wide), and second, by split- 
ting the sample within those bins based on morphology into 
early versus late types. Consequently, we use the small cor- 
rections to the best-fit ma sses determined from Table 1 in 
IMandelbaum et al.l ll2005lJ) for the case of no scatter in the 
M{L) relationship. Details of these corrections will be given 
in 



3 DATA 



The d ata used here are obtained from the SDSS llYork. et alJ 
I2OOOI) . an ongoing survey to image roughly tt steradi- 
ans of the sky, and follow up approxima tely one million 
of the detected objects spectroscopically ([Eise nstein et alJ 
I2OOII: iRichards et"alll2002l : IStrauss et al J 120021) . The imag- 
ing is carried ou t by drift-scanning the sky in ph oto- 
metric conditions llHogg et alJl200ll: llvezic et aljl20M) . in 



five bands (ugriz) iFukugita et al.lll996l: ISmith et al1l2002h 
using a specially-designed wide-field camera iGunn e*^!] 
Il998h . These imaging data are the source of the Large- 
Scale Structure (LSS) sample that we use in this pa- 
per. In addition, objects are targeted for spectroscopy us- 
ing these data iBlanton et alJ l2003lJ) and are observed 
with a 640 - fiber spectrograph on the same telescope 
iGunn et alj 1200^ . All of these data are processed by 
completely automated pipelines that detect and mea- 
sure photometric properti es of objects, and astrometri- 
cally cali brate the data I' Lupton et al.l 1200 ll : IPier et al] 
|2003; Tu cker et alj boOS'l. The SDSS is well underway, and 
has had five major data releases jStoughtoi^et al" 200! 
Abazaiian et alJ l2003l. l2004l. l2005l : iFinkbeiner'et al. .200. 
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3.1 Lenses 

The galaxies used as lenses are those targeted as the DR4 
MAIN spectros copic sample (478 3 dcg^), with calibration 
as described in iFinkbeiner et all (t2004j . Here we describe 
the quantities used for classifying the lenses into subsam- 
ples: morphology, stellar masses, luminosities, and density. 
We note that all samples described here are flux-limited, 
not volume-limited, with the flux limit nominally r < 
17.77 (Petrosian, extinction-corrected) but actually varying 
slightly in a known way across the survey area. Only those 
lens galaxies with redshifts 0.02 < z < 0.35 were used. 

An important feature of this work is the splitting of lens 
samples based on morphology within each stellar mass or lu- 
minosity subsample, since predictions of the quantities we 
are measuring may differ significantly for spiral and ellipti- 
cal galaxies. This split is carried out in practice by requir- 
ing that the parameter frac_deV output from the Photo 
pipeline be > 0.5 for early types, < 0.5 for late types. 
This parameter is determined by fitting the galaxy profile 
to frac_deV times the best-fit deVaucouleurs profile plus 
(1— frac_deV) times the best-fit exponential profile (and re- 
quiring <frac_deV< 1) in each band separately. For this 
work, to reduce the noise, we use the unweighted average 
over the g, r, and i bands. I n practice, fracdeV is highly 
correlated with galaxy colors. IStrateva et alJ (|20^) studied 
the effectiveness of morphological classification using either 
colors or de Vaucouleurs and exponential profile likelihoods 
for a sample of spectroscopic galaxies for which morphologi- 
cal classification was determined using spectral indices. The 
result is that for early types, selection by r-band frac_deV 
results in 96 per cent completeness, and is 76 per cent re- 
liable (that is, 96 per cent of the early types were indeed 
classified as early types, and 76 per cent of those classified 
as early types were actually early types). For late types, the 
completeness was 55 per cent, and the classification was 90 
per cent reliable. However, since we use the average of the 
results in g, r, and i, rather than just the results in r, and a 
more recent version of Photo, it is likely that the results of 
classification using frac_deV for this paper are even better 
than those quoted there. We note that selection by u — r 
color yielded 98 per cent completeness and 83 per cent re- 
liability for early types, and 72 per cent completeness and 
96 per cent reliability for late types, a slight improvement 
in the results over classification using profiles in that work. 
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though our use of profile information from three bands and a 
more recent version of Photo may decrease that advantage. 

Stellar mass estimates were obtained from the spectra 
using the same techniques as in [Kauffmann ct al. (200j), 
but for the full DR4 sample. A library o f 32,000 model star 
forma tion histories generated using the lBruzual fc Charlod 
l)2003h population synthes is models and the measured 
D„(4000) and H5a indices jTremonti et al.ll20o3 ) are used 
to obtain a median likelihood estimate of the 2-band mass- 
to-light ratio for each galaxy. By comparing the colour pre- 
dicted by the best-fit model to the observed colour of the 
galaxy, the attenuation of the starlight due to dust is also 
estimated. The corrected z-band magnitude and the MjL^ 
together yield an estimate of the stellar mass. These esti- 
mates have 95 percent confidence ranges of around ±40 per 
cent. Comparison of these measures of stellar mass, which 
are better measures of the mass in stars than luminosity, 
against the best-fit halo masses will allow us to learn about 
the efficiency of conversion of baryons in the halo to stars. 
We define this efficiency via the relation 

^ _ Mstellar MgteUar ^^-^ 

Mcent fib Mcentfb 

where fb is the cosmological baryon fraction. We note that 
this quantity is not the same as the traditional definition of 
star formation efficiency, the rate at which stars are forming 
compared to the available mass or surface mass density in 
HI or molecular gas. Hence, we use the term conversion ef- 
ficiency rather than star formation efficiency when describ- 
ing these results. There are two differences between these 
quantities. First, the star formation efficiency denotes the 
fraction of those baryons that have cooled off and formed 
part of the galax;y that then were included in stars. On the 
other hand, as shown in Eq.Q the conversion efficiency as- 
sumes that the baryon fraction in the galaxy is equal to the 
cosmological one, thus ignoring the fact that some signif- 
icant fraction of baryons actually are not in galaxies, but 
are located in, e.g., the hot gas in clusters, or warm-hot gas 
on the outskirts of galaxies. Furthermore, Eq. |7| includes all 
baryonic matter, such as helium, whereas the star forma- 
tion efficiency is defined relative to mass of hydrogen only. 
Because of these differences in definitions, our results for 
baryon conversion efficiency are also lower limits on the star 
formation efficiency. 

Table gives information about the stellar mass bins 
used for this work, each of which is a factor of two wide: the 
numbers of galaxies in each bin, the mean redshift within 
each bin (determined using the weights from the lensing 
analysis), and the fraction of galaxies classified as spirals 
{f spiral)- We note that because of the statistical error on 
each estimate, the bins are actually equivalent to tophats a 
factor of two wide in stellar mass convolved with approxi- 
mate Gaussians with a / (Msteiiar) ~ 0.2. Throughout this 
work, stellar masses Msteiiar are given in units of Mr?), hav- 
ing computed with h = 0.7 as in iKauffmann et all 

We also must consider systematic uncertainties in our 
results due to assumptions that are involved in the stel- 
lar mass determination . A full discussion can be found in 
IKauffmann et al.l J2003h : the main uncertainty that will af- 
fect our results is the uncertainty due to the IMF, since 
we must assume that the locally observed IMF applies uni- 



Table 1. The stellar mass subsamples used in this analysis, in- 
cluding mean weighted redshifts within each bin, the total number 
of galaxies, and the fraction of spirals. 



Sample 


^stellar 

[IOIOMq] 


^aai 




{^stellar) 

[IOIOMq] 


f spiral 


sml 


[0.5,1.0] 


23 474 


0.060 


0.74 


0.74 


sm2 


[1.0,2.0] 


40 952 


0.070 


1.5 


0.60 


sm3 


[2.0,4.0] 


66 503 


0.085 


2.9 


0.46 


sm4 


[4.0,8.0] 


90 019 


0.11 


5.7 


0.32 


sm5 


[8.0,16.0] 


82 734 


0.13 


11.0 


0.20 


sm6 


[16.0,32.0] 


39 729 


0.16 


21.0 


0.11 


sm7 


[32.0,64.0] 


8 096 


0.19 


40.0 


0.05 



versally. Changing from the iKroupal ||2Q^) IMF to another 
IMF (assuming that it is uniform across all samples), such 
as the Salpeter IMF, may rescale all the stellar mass val- 
ues, and therefore the derived conversion efficiencies, by a 
fixed value of up to 3 per cent. However, we note that 
ICappellari et all ||2QQ^ found with a sample of > L* ellip- 
ticals that the measured M/ L values appear inconsistent 
with the Salpeter IMF, which predicts too much mass, so 
this value of 30 per cent is likely conservative. To address 
the question of whether the IMF really is uniform across 
all samples, we note that lKroupa| i200ll) found no evidence 
for significant vari ations in the IMF within our own galaxy; 
assuming that the lCappellari et al.l (|22Q^ results represent 
typical elliptical galaxies, and the MW is a typical spiral, we 
can then infer that the IMF of spirals and ellipticals is not 
significantly different. Thus, even if the global IMF leads to 
some overall rescaling of the stellar mass values, many of 
the results which we will present are in the form of trends of 
conversion efficiency with stellar mass, or comparisons be- 
tween two different lens samples at fixed stellar mass, which 
are still valid regardless of this systematic uncertainty. 

Luminosities were determined using the r-band Pet- 
rosian apparent magnit udes (extinction-corr ected using the 
reddening maps from ISchlegel et alJ ll99Sh . k-corrections 
to z = 0.1 (the sam ple median) from kcorrect v1_11 
jBlanton et al.ll2003cl). p assive l uminosity evolution correc- 
tion frorn ^lanton^tall i2003lJ) . and the distance modulus 
determined using ft = 1, yielding Mr = m + 1.6{z — 0.1) — 
{K + DM). These luminosities m ay be underestim ated at 
typically the 10-20 per cent leve l llWest et al.ll2005|) due to 
the sky subtraction systematic iAdelman-McCarthv et all 
120051: iMandelbaum et a^]|2005a^ . an overestimate in the lo- 
cal sky estimate within ~ 90" of bright (r < 18) galaxies and 
stars, thus affecting the entire lens sample at some level. The 
actual underestimation varies with apparent magnitude and 
size, and changes the estimated M/L; since the effect has 
not been well-quantified for the full sample and is within 
the 95 per cent CL, we have applied no correction. 

One possible systematic in the comparisons between 
stellar mass and luminosity is the fact that the luminosity 
is the full Petrosian magnitude (which contains essentially 
all the light of an exponential profile, and about 80 per cent 
of the light in a de Vaucouleurs profile), whereas the stellar 
mass is estimated using spec tra measured using a 3" aper- 
ture. |Kau^^^n^^^L|j200^ discuss this "aperture bias" in 
detail, and while it is found to infiuence the results at some 
level, the variation in Msteiiar /L due to this systematic are 
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Table 2. The luminosity subsamples used in this analysis, includ- 
ing mean weighted redshifts and luminosities within each bin, the 
total number of galaxies, and the fraction of spirals. 



Sample 


Mr 








(L/L,) 


f spiral 


LI 


-17 > Mr > 


-18 


10 047 


0.032 


0.075 


0.80 


L2 


-18 > Mr > 


-19 


29 730 


0.047 


0.19 


0.69 


L3 


-19 > Mr > 


-20 


85 766 


0.071 


0.46 


0.53 


L4 


-20 > Mr > 


-21 


141 976 


0.10 


1.1 


0.35 


L5f 


-21 > Mr > - 


21.5 


60 994 


0.14 


2.1 


0.23 


L5b 


-21.5 > Mr > 


-22 


34 920 


0.17 


3.2 


0.16 


L6f 


-22 > Mr > - 


22.5 


13 067 


0.20 


4.9 


0.08 


L6b 


-22.5 > Mr > 


-23 


2 933 


0.22 


7.7 


0.05 



insignificant compared to its variation with galaxy proper- 
ties such as luminosity and profile type (parametrized there 
using the concentration index). 

In order to obtain mass-to-Ught ratio s, solar lumi- 
nosities were determined using results from iBlanton et alJ 
l|2nn3lJ^ with M, = -20.44 and Msoiar = 4.76, yielding 
luminosities in /i~^Lq (L* = 1.2h-^ X 10^°Lq). Relevant 
information for luminosity subsamples is shown in Table |5| 
including numbers of galaxies, mean weighted redshifts and 
luminosities, and spiral fractions. We note that we split 
the two brightest bins into half-magnitude bins, because 
we would like the opportunity to better constrain the vari- 
ation of mass with light in that regime, and because the 
central halo mass distribution may become wider at higher 
luminosities since some of the galaxies are Brightest Clus- 
ter Galaxies (BCGs) of clusters with large halo masses, and 
others are field galaxies hosted by smaller halos, giving a 
broader halo mass distribution at the high- luminosity end. 

Finally, we need a measure of the local galaxy envi- 
ronment. Many estimators, including nth nearest neighbor 
(3-d or in projection), counts in an aperture (again, 3-d or 



literature (for example, iRamella et al. 


I2OOII: 


Marinoni et al 


2002'; 'Hoee et al."'2003'; Baloah et al.| 


2Q()4b 


: iBlanton et al. 


20I)5b.; ,CaQner_et al. 2005). Here, we choose a verv simple 



one, spectroscopic galaxy counts in cylinders of radius 1 
h'^Mpc and line-of-sight length Av = ±1200 km s"\ These 
numbers are compared to the number of random points in 
the same cylinders, thus taking into account the angular and 
radial variations in number density, survey boundaries, and 
other issues. While galaxies excluded in the spectroscopic 
survey because of fiber collisions are not used for the stel- 
lar mass study (since they lack spectra), they are used for 
the density estimates in order to avoid underestimating the 
environment measure in rich clusters. The galaxies without 
spectra due to fiber collisions are given redshifts equal to 
those of the nearest neighbor. We note that because the en- 
vironment measurement requires a careful knowledge of the 
survey completeness as a function of position, only those 
galaxies included in the area covered by the LSS DR4 sam- 
ple from the NYU V alue-Added Galaxy Catalog (VAGC; 
IBlanton et"ai]|2005al) were used to obtain environment esti- 
mates. 

Our estimates were derived using 20 times as many ran- 
dom points as real galaxies; we note that for higher redshifts, 
these environment estimates can be quite noisy or may fail 
entirely, because the radial selection function leads to a low 



density of objects. Furthermore, galaxies at the lower red- 
shift limit or near survey boundaries may not have envi- 
ronment estimates if no random galaxies were found in the 
cylinders around them; thus, when splitting the sample at 
the median density within each stellar mass or luminosity 
bin, only some fraction of the sample was used (ranging from 
93 per cent for bins at L <~ L*, down to 40 per cent for the 
brightest galaxies) . Fortunately, the lensing signal itself pro- 
vides a reasonable test of the environment estimate, because 
we can see whether the signal at 1-2 /i^^Mpc scales is con- 
sistent with the lens sample primarily being in the field or 
in groups and clusters. At fixed luminosity or stellar mass, 
the spiral sample (defined by frac_deV < 0.5) had a lower 
median and mean environment estimate than ellipticals, as 
expected from previous studies cited in ^ 

Because of the need for spectra to determine stellar 
masses, our sample must exclude galaxies for which spec- 
tra were not obtained due to fiber collisions. Fibers cannot 
be placed closer than 55" (~ 80 /i~^kpc at z ~ 0.1), so if two 
targets are closer than this separation, only one will have a 
spectrum. This restriction is alleviated in roughly one third 
of the sky by the use of overlapping plates, but nonetheless, 
roughly 7 per cent of targets do not have spectra. This loss of 
targets is naturally worse in high-density regions, and will 
therefore tend to affect the satellite contribution, decreas- 
ing it in a scale-dependent way, and therefore changing its 
shape. However, as our results will show, we are not highly 
sensitive to the shape of the satellite contribution. 



3.2 Sources 

The source sample is the same as that from 
iMandelbaum et all l2005al (hereinafter M05), who ob- 
tained shapes for more than 30 million galaxies in the 
SDSS imaging data down to magnitude r = 21.8 (i.e. four 
magnitudes fainter than the SDSS spectroscopic limit). 
This section briefly describes the M05 pipeline, also known 
as Reglens. 

The M05 pipeline obtains galaxy ima ges in the r and 
i fllt ers from the SDSS "atlas images" dStouehton et all 
l2002h . The basic principle of shear measurement using these 
images is to fit a Gaussian profile with elliptical isophotes 
to the image, and define the components of the ellipticity 



{e+,ex) 



1 



{b/o 



1 + (b/a) 



-(cos 2<jf>, sin 2(j}), 



(8) 



where b/a is the axis ratio and is the position angle of the 
major axis. This is then an estimator for the shear. 



(7+,7x) = — ((e+,ex)), 



(9) 



where TZ 0.87 is called the "shear responsivity" and repre- 
sents the response o f the ellipticity (Eq . O to a small shear 
jKaiser et alJ Il995l : iBernstein fc JarvislboO^I . In practice, 
a number of corrections need to be applied to obtain the 
ellipticity. The most important of these is the correction 
for the smearing and circularization of the galactic images 
by the PSF; M05 uses the PSF maps obtained from stel- 
lar images by the psp pipeline jLupton et alJ 1200111 . and 
corrects for these usin g the re-Gaussianization technique of 
iHirata fc SeliakI i2003l) . Re-Gaussianization corrects for the 
PSF while taking into account the non-Gaussianity both of 
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the PSF and of the galaxy profile. A smaller correction is for 
the optical distortions in the telescope: ideally the mapping 
from the sky to the CCD is shape-preserving (conformal), 
but in reality this is not the case, resulting in a nonzero 
"camera shear." In the SDSS, this is a small effect (of order 
0.1 per cent) which c an be identified a nd removed using the 
astrometric solution l|piereiai]|2003). Finally, a variety of 
systematics tests are necessary to determine that the shear 
responsivity TZ has in fact been determined correctly. We 
refer the interested reader to M05 for the details of these 
corrections and tests. 

M05 includes a lengthy discussion of shear calibration 
biases in this catalog; we will only summarize these issues 
briefly here. Our source sample is divided into three subsam- 
ples: r < 21 , r > 21, and high-reds hift Luminous Red Galax- 
ies (LRGs, lEisenstein" et al."200lD, defined using color and 
magnitude cuts as in M 05 using selecti on criteria related to 
those from lEisenstein et all feoOl ) and'Padmanabhan ct al.' 
l|2005() . Using simulations from iHirata fc Seliak (.2003.) to 
estimate the PSF dilution correction and analytical models 
for selection biases and other issues that affect shear cali- 
bration, we place 2a limits on the multiplicative shear cali- 
bration bias of [—5, 12] per cent for r < 21, [—8, 18] per cent 
for r > 21, and [—6, 19] per cent for LRGs. 

As shown in Eq. |21 the lensing signal AE is a product 
of the shear and factors involving lens and source redshifts. 
Since the lenses have spectroscopic redshifts, the primary 
difficulty is determining the source redshift distribution. We 
take three approaches, all described in detail in M05. For 
the r < 21 sources, we use photometric redshifts and their 
error distributions determined using a sample of galaxies 
in th e Groth strip with redshifts from DEEP2 LDa^^^^ 
l2003l: iMadgwick erai] 120031: ICoil et alJ l2004l : 1 Davis et al 
|200J), and require Zb > zi + 0.1 to avoid contamination 
from physically-associated lens-source pairs. For the r > 21 
sources, we use redshift distributions from DEEP2. For the 
high-redshift LRGs, we use photometric redshifts and their 
error distributions determined using data from the 2dF- 
Sloan LRG and Quasar Survey (2SLAQ), and presented in 
[Padmanabhan ct al. (2005). 

Finally, we have placed constraints on other issues af- 
fecting the calibration of the lensing signal, such as the 
sky subtraction problem, intrinsic alignments, magnification 
bias, star-galaxy separation, and seeing-dependent system- 
atics. As shown in M05 the calibration of the signal using 
the three source samples agrees to within 10 per cent, with 
a total la calibration uncertainty estimated at 7 per cent 
(r < 21) or 10 per cent (r > 21 and LRG). 



3.3 Signal computation 

Here we describe the computation of the lensing signal. Lens- 
source pairs are assigned weights according to the error on 
the shape measurement via 



lensing signal in 46 logarithmic radial bins from 20 h ^kpc 
to 2 /i~^Mpc as a summation over lens-source pairs via: 



Wis 



aj + at 



(10) 



where cr|jv, the intrinsic shape noise, was determined as a 
function of magnitude in M05, figure 3. The factor of 
downweights pairs that are close in redshift. 

Once we have computed these weights, we compute the 



AE(i?) 



Wis 



(11) 



where the factor of 2 arises due to our definition of ellipticity. 

There are several additional procedures that must be 
done when computing the signal (for more detail, see M05) . 
First, the signal computed around random points must be 
subtracted from the signal around real lenses to eliminate 
contributions from systematic shear. Second, the signal must 
be boosted, i.e. multiplied by B{r) = n{r) / n,.and{r) , the ra- 
tio of the number density of sources relative to the number 
around random points, in order to account for dilution by 
sources that are physically associated with lenses, and there- 
fore not lensed. 

In order to determine errors on the lensing signal, we di- 
vide the survey area into 200 bootstrap subregions, and gen- 
erate 2500 bootstrap-resampled datasets. To determine er- 
rors on fit parameters, we perform the fits on each bootstrap- 
resampled dataset, and use the distribution in parameter 
space to determine confidence intervals. 



4 RESULTS 

4.1 Lensing signal 

In Fig. we show the lensing signal in stellar mass bins 
with la errors, for spiral (frac_deV< 0.5) and elliptical 
(f rac_deV> 0.5) galaxies, with the best-fit halo model signal 
(with parameters to be described in H4.2II . Fig. |5| shows the 
same for luminosity bins; the signal was noisy enough for 
the brightest bin (L6b) that errors could not be determined 
on fits, so that bin is not shown. Finally, Fig. |3 shows the 
results in luminosity bins for low- and high-density samples 
of ellipticals only. The results have been rebinned for eas- 
ier viewing; correlations between radial bins are minimal at 
small radius, and reach a maximum of about 0.2-0.3 at the 
outermost bins for the lower luminosity or stellar mass bins 
(at lower z) and are negligible at all radii for the brighter 
bins. While the model fits are shown on the plots, we defer 
discussion of the best-fit parameters to the following sec- 
tions, and focus here on a comparison of the lensing signal. 

First, we consider Fig. There appears to be some 
detection of signal on all scales for all but the lowest bin, 

5 X 10^ < MstMar/MQ < 10^^ For sm2-sm4, comprismg 
(1 — 8) X 10^" Mq, the signal for early and late types appears 
statistically consistent for r <~ 100 /i~^kpc, indicating that 
the average halo masses (and therefore conversion efficien- 
cies) for the two types of galaxies are similar. For sm5 and 
sm6, comprising (8 — 32) x 10^'' Mq, the signal on small 
scales is larger for early types than for late types, indicating 
a larger average halo mass for early types, so stellar mass is 
no longer a good tracer of halo mass for this range of stellar 
masses. This result is not entirely surprising, since for these 
bins, some fraction of the red galaxies are LRGs that may be 
the BCGs of clusters, so the halo masses determined using 
lensing are actually the masses of the entire cluster rather 
than of individual galaxies. On larger scales, the signal for 
early types is higher than for late types for each stellar mass 
bin to a degree that is highly statistically significant, consis- 
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AS for stellar mass bins 




100 1000 
r [kpc/h] 



Figure 1. AS in stellar mass bins for early (red hexagons, solid line) and late-type (blue triangles, dashed line) galaxies. For the highest 
stellar mass bin, only the signal for early-types is shown since they are 95 per cent of the sample. All errors are la. 



tent with previous findings that red galaxies reside in denser 
environments than blue galaxies. 

Next, we consider Fig. |21 The detection is marginal in 
LI and L2, comprising magnitudes — 19 < Mr < —17. For 
the bins with a robust detection, L3 (with (L/L^) = 0.46) 
is the only one for which the signal for early and late types 
appears to be highly consistent at small scales; at higher 
luminosity, it appears to be higher for early types. On large 
scales, the signal for early types is again larger than for late 
types for all luminosity bins, indicating a tendency for early 
types to be in denser environments. 

Finally, we consider Fig. |21 which shows AE in lumi- 
nosity bins for early-type galaxies only, with the sample di- 
vided at the median environment measure in that bin. We 
remind the reader that environment estimates were not ob- 



tained for a large fraction of the galaxies in the brightest 
bins due to the limitations of the environment estimator, so 
few lenses were used for those bins. For L1-L5, it is appar- 
ent from the signal on r > 300 /i^^kpc that the environment 
split efficiently separates the galaxies into those that reside 
in overdense versus underdense regions. This separation is 
particularly marked for L3 and L4, with the signal on 1-2 
/i~^Mpc scales differing between the high- and low-density 
samples by more than an order of magnitude. For the bright- 
est bins, the difference is less striking; it is possible that our 
environment estimator is simply too noisy at higher red- 
shifts, or that enough of the L6 early-type galaxies reside in 
groups and clusters that the environment simply does not 
vary much across the bin, leading to a small variation in 
signal. For all bins, the signal on small scales does not vary 



AE for luminosity bins 
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Figure 2. AS with best-fit lialo model in luminosity bins for early (red hexagons, solid line) and late-type (blue triangles, dashed line) 
galaxies. All errors are Icr. 



significantly with the environment. As we show below, the 
L* high density sample is dominated by satellites. We see 
that the tidal stripping of dark matter around satellites in- 
side groups and clusters cannot be maximally efficient, since 
this would have been seen as a suppression of signal on small 
scales. 



4.2 Halo model fits 

Here we discuss the goodness-of-fit for the halo model fits 
before proceeding to discussions of the best-fit parameters. 
However, we note that the values of are not expected 
to follow the usual distribution b ecause of the noisin ess 
of bootstrap covariance matrices (see iHirata et alJl20o3 for 



a fuller description of these results). For our fits with 40 
degrees of freedom and 200 bootstrap regions, the expected 
value of is 50.4, not 40. It is clear from the figures that 
the halo model signal generally is a close match to the data; 
for example, for early types, the for the fits for the seven 
stellar mass bins were 37.5, 39.9, 31.5, 36.5, 46.9, 40.5, and 
42.0 respectively, with p- values (i.e. the probability to exceed 
this value by chance if the model is a realistic description 
of the data) of 0.85, 0.79, 0.96, 0.87, 0.58, 0.78, 0.73. The 
results may imply that the bootstrap errorbars are slight 
overestimates. Results for late types in stellar mass bins, 
and the splits by luminosity, were similar. 

One concern regarding these fits is that for the highest 
stellar mass and luminosity bins the halo model underes- 
timates the signal by a significant amount on small scales 
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AS for ellipticals in luminosity bins (split by environment) 
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Figure 3. AS with best-fit iialo model in luminosity bins for early-type galaxies, divided into overdense (magenta hexagons, dashed 
line) and underdense (black triangles, solid line) samples. All errors are Itr. 



(r < ~ 40 ^kpc). Because of the sky subtraction prob- 
lem Adelman-McCarthv et al.l 120051: iMandelbaum et alJ 
l2005air only scales r > 30 /i ^kpc are used for the fits even 
though we have plotted down to 20 /i^^kpc, so this discrep- 
ancy between the best-fit and observed signal on these scales 
does not cause a major increase in the x^; however, it is still 
a concern, since the magnitude of the effect is larger than 
the estimated effect due to sky subtraction (at most 10 per 
cent). Because the large magnitude of this discrepancy does 
not seem attributable to any of the small-scale systematics 
such as intrinsic alignments, the sky subtraction problem, 
or magnification bias, the discrepancy seems to suggest a 
problem with the assumed profile itself, likely due to the 
efl'ects of stellar component and the associated dark matter 
contraction caused by it, which we have completely ignored. 



A detailed study of this discrepancy between the model and 
the observed signal will be investigated in the context of 
dark matter profile constraints from g-g lensing in future 
work. 

O ne final issue we must consider is that brought up in 
IMand elbaum ct al. ( 2005h) of the corrections to the central 
halo masses due to broadness of the central halo mass distri- 
bution, which naturally has some width due to the width of 
the stellar mass or luminosity bins, and which may be fur- 
ther broadened by scatter in the mass-luminosity (or mass- 
stellar mass) relationship. As noted there, the best-fit central 
halo mass underestimates the mean value (and overestimates 
the median, which is not of interest in this work) by some 
amount due to this width. In this work, we use the correc- 
tions in the "no-scatter case" which are purely due to the 
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width of the bins that lead to a wide central halo mass dis- 
tribution, because we anticipate that our division into early 
and late type samples will lead to narrower central halo mass 
distributions. Those corrections, as described there, are to 
increase the best-fit masses by 6 per cent, 11 per cent, and 
36 per cent for L3, L4, and L5 (1 magnitude wide), respec- 
tively. We extrapolate these values to higher luminosity and 
use no correction for L1-L2, 6 per cent, 11 per cent, 25 per 
cent, 50 per cent, and 66 per cent for L3, L4, L5f, L5b, and 
L6f, respectively. Fundamentally these values are somewhat 
uncertain, and indeed represent one of the main modeling 
uncertainties of this paper. All corrections have been applied 
to values in the tables and text in the following sections. For 
stellar mass samples, we use the average luminosity of each 
bin to find a correspondence between stellar mass and lumi- 
nosity bins, and use corrections of 0, 6, 8, 11, 22, 45, and 66 
per cent for stellar mass bins 1-7, respectively. While these 
corrections introduce some uncertainty in the masses them- 
selves, they do not affect comparisons between morphology 
samples at the same stellar mass or luminosity since the 
same corrections are applied to each morphology sample. 

4.3 Central halo masses 

Here we present results derived from the best-fit central halo 
masses and their 95 per cent confidence intervals. Due to 
the imposition of the constraint Mcent > 0, the error distri- 
butions at low luminosity and stellar mass are highly non- 
Gaussian, and the confidence intervals are determined using 
the distribution of Mcent values for the bootstrap subsam- 
ples. 

First, we consider the relationship between halo mass 
and stellar mass. Table|H]shows the best-fit halo mass and its 
95 per cent confidence interval as a function of stellar mass 
for early and late type galaxies. For reference, the average 
r-band luminosity for each sample is also shown. Figure |3] 
shows a plot of halo mass as a function of stellar mass, and 
of central galaxy conversion efficiency as a function of stellar 
mass, where we define conversion efficiency -q via Eq. Qwith 
Qm = 0.27, = 0.046, and h = 0.7. 

Table O and Fig. 2] reveal several interesting trends 
in the relationship between stellar mass and halo mass. 
First, stellar mass is a good tracer of halo mass for 
M stellar < 10^^^ Mq (corresponding roughly to L < 1.5L*); 
that is, both early- and late-types have the same average 
relationship between halo and stellar mass in this regime. 
Above this value of stellar mass, the central halo mass for 
early types is larger than that for late types. One possible 
explanation for this trend is that early types are far more 
likely to be the central galaxy of a group- or cluster-mass 
halo than late types, which are only at the center of smaller 
halos. Thus, the Mcent values in this table and figure for 
early types at higher values of stellar mass are likely to re- 
flect the mass of the entire group/cluster rather than that 
of the central galaxy only. 

We also consider the lower panel of Fig. 0] the central 
galaxy conversion efficiency 77 as a function of M stellar- For 
those bins with a marginal detection of signal, we are able 
to place only weak constraints on 77, which may at times 
even have a central value outside of the range of reason- 
able values < r; < 1; however, in all cases the 2a errors 
shown in the table and figure include part of the range of 



Table 3. The central halo mass determined for stellar mass sub- 
samples separately for early and late type galaxies, and the re- 
sulting conversion efficiencies determined according to Eq. |7| All 
confidence intervals are 95 per cent. 
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Figure 4. Central halo mass (top) and conversion efficiency in 
central galaxy 77 (bottom) as a function of stellar mass for early- 
type (red hexagons, solid line) and late-type (blue triangles, dot- 
ted line) galaxies. All errorbars shown are 95 per cent confidence 
interval. 



reasonable values. The particularly useful part of this exer- 
cise is the lower-limit in r) derived from the upper-limit on 
the halo masses. For Msteiiar > 10^^^ Mq, we see a reflection 
of the trend noted for the early types in the upper panel, 
that many of them are central galaxies of groups or clus- 
ters, and consequently the interpretation of Msteiiar / Mcent 
in terms of conversion efficiency is no longer reasonable. (To 
do so, we would have to include the stellar masses of all other 
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Figure 5. Central halo mass (top) and Meant/ L as a func- 
tion of luminosity for early- type (red hexagons, solid line) and 
late-type (blue triangles, dotted line) galaxies. Dashed lines also 
show the best-fi t powe r-law M{L) for early and late types from 
iGuzik fc Sell^ |2 



as discussed in the text. All errorbars 



shown are 95 per cent confidence intervals. 



group/cluster members or include the mass of the central 
galaxy only rather than the full cluster mass.) However, in- 
cluding stellar mass of other group/cluster members would 
tend to raise rj, so while the central value of rj is likely to 
be underestimated at high Msteiiar, the lower limits given 
here are robust. Hence, for late type galaxies, we can con- 
clude that for typical galaxies with MstMar ~ 10^^ Mq, the 
conversion efficiency is typically 35 per cent, with a lower 
limit of 18 per cent (95 per cent CL). For early types, the 
typical galaxy has a conversion efficiency of only 14 per cent, 
or a lower limit of 11 per cent. These results imply a differ- 
ence in conversion efficiency between typical early and late 
type galaxies of roughly a factor of two or more, with the 
efficiency for both early and late types declining at higher 
stellar masses. 

Next, we consider the relationship between r-band lu- 
minosity and halo mass. Table|l]and Fig.|S]give the relevant 
parameters. This table and figure reveal several interesting 
points about the relationship between r-band luminosity and 
halo mass. First, the luminosity is a good tracer of central 
halo mass for L <~ 0.8Z/«, with early- and late- types below 
this luminosity having consistent halo masses. At higher lu- 
minosities, we see a similar trend as for the split by stellar 
masses, with larger central halo masses at the bright end for 
early types, possibly explained by them being at the centers 
of group- or cluster-mass halos. Hence, while the relation- 
ship between mass and luminosity appears to be consistent 
with a single power-law for late types, the same is not true 
for early types, due to the steepening in this relationship 
at the bright end (which, in this regime, is more accurately 
a relationship between the BCG luminosity and the clus- 
ter mass rather than an individual galaxy luminosity and 
mass). In the lower panel of Fig. |S1 we note that because 
of this inclusion of the mass in all group/cluster members. 



Table 5. The central halo mass determined for luminosity sub- 
samples separately for early-type galaxies in low- and high-density 
regions. All confidence intervals are 95 per cent. 



(L/L.) M,e„t [IQH/i-IMq] 
Low density 



High density 



0.075 

0.19 

0.47 

1.1 
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4.9 
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12.0 
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-1.9 
+8 



+30 
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Figure 6. Central halo mass as a function of luminosity for early- 
typo galaxies in low-density (magenta hexagons, solid line) and 
high-density (black triangles, dotted line) samples. All errorbars 
shown are 95 per cent confidence intervals. 



the central values and lower limits of Mc^nt/L are not as 
relevant as the upper limits. Consistent with the results de- 
rived using stellar masses, Mcent/L is larger for early types 
than for late types by about a factor of two around L* , and 
the value of M^ent/L decreases significantly for higher lu- 
minosities. As for the relationship between halo mass and 
stellar mass for isolated halos, we need to use environment 
estimates to explore what happens to field galaxies at high 
luminosities in order for the values for early types to have a 
simple physical interpretation. 

Next, we consider the relationship between local density 
and halo mass for ellipticals, where we divide each luminos- 
ity bin at the median density in that bin. Table|K|and Fig.|S| 
contain the results. 

The signal in Fig. |H] and the best-fit masses in Ta- 
ble |K| and Fig. |S| appear to suggest that early-type galaxies 
in both low-density and high-density regions have consis- 
tent central halo masses, particularly for L <~ L* where 
the high-density sample is exclusively satellites, rather than 
BCGs. Another interesting point to note is that discussed 
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Table 4. The central halo mass determined for luminosity subsamples separately for early and late type galaxies. All confidence intervals 
are 95 per cent. 



Mr 


(L/L,) Meant Meant /L 
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-22) 


5.0 


O'o_174 
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"'^-203 
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in iMandelbaum et al.l ll2005bft . that a large amount of tidal 
stripping of satellite profiles in clusters would cause a clear 
depression in the lensing signal on scales below the virial ra- 
dius. For example, a scenario predicted from a combination 
of N-body simulations w ith semi-analytic m odels of galaxy 
formation suggested by iGao et all yOM^) is that a very 
large fraction may have nearly all the dark matter stripped. 
(As stated in ^ we do assume that all satellites have 50 per 
cent of the matter stripped.) No such major depression in 
the signal relative to that for field galaxies is seen on 50-100 
/i~^kpc scales, and therefore we can rule out scenarios that 
have most satellites strongly stripped. 



4.4 Comparison against past studies 

Other works have attempted to determine whether halos of 
satellite galaxies arc tidally stripped. The methodology pre- 
sented by .Nataraia n & Kncib ( 1997) allowed for the mass 
distributions of individual cluster galaxies to be mapped out, 
in addition to the large- scale cluster m ass distribution. This 
method was applied bv lNataraian et al. (1998), who found 
some suggestion that tidal stripping may affect E galaxies 
less than 5*0 galaxies (based on diff erences in velocity dis- 
persions), and lNataraian et al.||2002| (with a lo nger descrip- 
tion of methodoloev in lNe^iraian et al.ll2004) . who found 
clear evidence for tidal stripping, with field-scale halos for 
Lf early types excluded at the lOcr level. Their density pro- 
file for the galaxy-scale halos differed from the one used in 
this work, but shares the property that below the truncation 
radius, AE oc (our data is consistent with a power- law 
index of —0.9, so quite similar), and oc beyond trun- 
cation radius rt, which is the same as our model, but with 
the transition happening gradually rather than abruptly as 
in our model. 

For a typical cluster, the truncation radius for L, 
early type galaxies was found in Nataraian et al. (2002) to 
be 40 kpc, or 28 fe-^kp c; our best-fit value of Mcent = 
151* X W^^ Mq for high-density early type galaxies at 
L/Lt — 1.1 indicates a virial radius of r = 2 88 /i~^kpc, 
which combined with the lNataraian et alj i2002 t) result sug- 
gests no dark matter beyond 10 per cent of the virial radius 
for an equivalent halo in the field, or 85 per cent of the halo 
stripped (rather than our assumption of no matter beyond 
40 percent ofrvir, or half the mass stripped). Tidal stripping 
beyond 28 /i~^kpc would give a much more suppressed signal 



on scales ~ 28 to 200 ft^^kpc (above which the group/cluster 
contribution dominates) than what we see here. However, 
we note that those results are for a small number of mas- 
sive clusters (five), whereas our results are an average over 
all group- and cluster-sized halos. Furthermore, our results 
average over all separations from the cl uster center, whereas 
the results from lNataraian are all within 50-100 

arcsec, or a characteristic scale of 300 /i~^kpc, significantly 
smaller than the virial radius of a massive cluster, and there- 
fore in the regime where the tidal stripping is expected to 
be more efficient, so it is not possible to draw quantitative 
conclusions from this comparison. Gavazzi et al. ( 200^ also 
found, using a weak lensing analysis of a single superclus- 
ter, that for early types, the typical scale of the dark matter 
halos of satellites decreased from about 300 kpc on the out- 
skirts of clusters to 200 kpc near the cluster cores; this result 
is consistent with our assumed model of tidal stripping, and 
with our findings. 

A number of previous studies have attempted to 
determine the relat ionship between mass and light. 
iGuzik fc SeliakI (|20o2) found a roughly power-law M(L) = 
Mt{L/Lt)^ in each band using a halo model analysis of 
g-g weak lensing data in each band from the SDSS, but 
with a much smaller galaxy sample, and different meth- 
ods of calibrating the shear and the redshift distributions 
than are used here. Their results in the r-band were M = 
(8.96±1.59) X 10" (L/L*)^'^^±° °''/i-^ Afg, as determined us- 
ing luminosity samples spanning the range 0.6 < L/L* < 7. 
However, these results assumed a = 0; using a more reason- 
able value of a ~ 0.2 gives values of M* that are 10-20 per 
cent lower. Furthermore, that work defined the mass using 
200perit rather than 180p, which means that in order to com- 
pare their M« with ours, we must increase their value of M* 
by approximately 30 per cent. Thus, from the combination of 
these two effects, we increase their Af* by 10 per cent, yield- 
ing Af. = (9.9±1.7) X lO^^fe"^A^0 (Icr error). In TableU we 
can see that for the luminosity bin closest to L* , our central 
halo masses are (14.9l4;6) and (T.lla;*) x lO"/i"^Af0 for 
early- and late-types respectively. So, using the 35 per cent 
spiral fraction for this bin from Table |21 we get an average 
halo mass of (12.212:5) x 10"/i"^AfQ {2a), which is 20 per 
cent higher than but statis tically consistent with the value 
from lGuzik fc SeliakI (|2002h . We also compare against their 
M, values for early- and late-types in the r-band, also cor- 
rected upwards by 10 per cent, which are (11.8 ± 2.8) and 
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(3.6 ± 2.3) X 10"/i"^Mq (where the fits allowed a to vary, 
unlike the above estimates, and fixed /3 to 1.51), which is 
also lower than but consistent with our results. In Fig. |^ 
the top panel includ es as dashed lines the M{L) for early 
and late types from iGuzik fc Se liak (2002); as shown, the 
fit, while not perfect, is not totally ruled out by our 2a er- 
rors. In particular, for late types, the slope appears reason- 
able but the amplitude slightly low compared to our results, 
whereas for early types, the deviation between our results 
and the fits is most prominent for high luminosities, where 
the fits underestimate Mcsnt (i.e., our results are consistent 
with a steeper slope /3). We note that there are a number 
of other differences in the modeling between this paper and 
that one (e.g., the determination of a and a lack of correc- 
tions for width of the central halo mass distribution in the 
older work), and differences in calibration of the lensing sig- 
nal, so a more detailed comparison is not necessarily useful 
or meaningful. 

Next, we compare against the recent work by 
iHoekstra et alJ (|^05). A number of factors make direct 
comparison difhcult: the different average lens redshift (with 
complications due to luminosity evolution) , the use of differ- 
ent passbands (V^, B, and R rather than SDSS bands), the 
different definitions of virial mass, and the selection of only 
isolated galaxies to avoid the need to fit for a non-central 
term. However, a number of their findings may be compared 
directly against ours. First, they find M « L^'^ in all three 
bands, similar to lGuzik fc SeliaM i2002l) : the previous para- 
graph includes a discussion of this result relative to ours 
as shown in Fig. |3 Next, when accounting for luminosity 
evolutio n, difference in passban ds, and different definitions 
of mass, iHoekstra et al.l (|200^ find that t heir masses are 
about 25 per cent higher than those from iGuzik fc SeliakI 
i2002^ when comparing against the g-band results, but still 
statistically consistent. This result implies that our results 
for the masses are consistent with those from lHoekstra et alJ 
i200a|). since our masses ar e also slig htly higher than those 
from lGuzik fc SeliakI (|20o3). Finally, IHoekstra et"ai] (120051) 
find conversion efficiencies a factor of two higher for late 
types than for early types, 33 per cent versus 14 per cent. 
This difference in morphology appears to be consistent with 
our results (see Fig.llJ for stellar masses larger than IO^^AIq 
or L > I/,; for lower stellar masses or luminosities, we lack 
statistical precision to make concrete statements. 

The conditional lumi n osity function fits to 2dF data 
from Ivan den Bosch e t alj l|2 003l) vield M(L) that is quite 
similar to that given in lVal^^Qstrikeil i2004l) based on em- 
pirical models, which is M oc L°'^^ at low mass (or luminos- 
ity), and M (X L'^'® at the high mass end, where the luminos- 
ity is that of the BCG alone (i.e. not all the cluster galaxies 
combined) and the mass is that of the full cluster. While we 
are unable to constrain the power-law slope very well at low 
luminosities because of the errors, the three brightest lumi- 
nosity bins with (L5f, L5b, L6f) give M oc j;,2-7±0-6 (95 p 
cent C L), with exponent a bout 3a below t helVale fc Ostrikeil 
i2004l) model predictions. Ivan den Bosch et al.ril20o3n igure 
1 also shows that the conversion efficiency is highest for L* 
galaxies, with rapidly increasing M/L for lower and higher 

meisses, consistent with our results. 

Analytical models of the Milky Way iKlvpin et alJ 

|20o3) predict halo masses of 7 X lO"/i"^M0, where the 
mass is defined as that within the radius within which the 



average density is equal to 340p. To compare against our re- 
sults, their mass estimate must be increased by 15 per cent, 
giving 8 X 10^^fe~^MQ, which we compare against our result 
for L, late-type galaxies, Mcent = {7 ■'i-tl'^) x 10"/i"^Mq 
(95 per cent CL). Alternatively, we note that their table 2 
suggests a total stellar mass of 6 x 10^° Mq, which is typical 
of L, galaxies and which (according to our Table I^J gives a 
halo mass of Mcent = 14ly x 10^^h~^MQ. Thus, it appears 
that this analytical model of the Milky Way is consistent 
with our results at the 95 per cent CL. 

It is also worthwhile to compare against the lensing 
signal from A'^-body simulati ons. In this case, we com - 
pare against the fit results from lMandelbaum et alj ll2005lJ) . 
which showed the lensing sig nal for three luminosit y bins 
from simulations described in lTasitsiomi et al] (|2Q23)i both 
with and without scatter in the mass-luminosity relation- 
ship: L3, L4, and L5 (a bin one magnitude wide that in- 
cludes L5f and L5b). Without scatter, the values of Mcent 
(mean) in each bin were 5.1, 18, and 132 x 10^^/i~^Mq; 
with scatter, they were 7.5, 29, and 117 x 10"/i"^MQ. For 
this paper, if we combine the results in each bin (averaging 
the results with different morphologies using the spiral frac- 
tion for weighting purposes), we get 4.3 ± 2.5, 12 ± 3, and 
(82± 15) X lO"fe"^M0. Therefore, the results from simula- 
tions seem to give somewhat higher masses than are found in 
the real data, though just within the 95 per cent confidence 
intervals, except for L5. 

Finally, we comp are against the baryonic mass func- 
tion determination by iRead fc Trenthan] i2005ft that uses 
a variety of data sources. In that work, the fraction of 
baryons in galaxies is estimated to be ~ 10 per cent, and 
of those, ~ 80 per cent are in stars (with the star forma- 
tion efficiency varying b y morphological type). Hence, the 
iRead fc TrenthamI ^2005^ results suggest average rj values of 
0.08. While our peak values of rj (around W^^ Mq) are 
higher than this, the abundance of galaxies at much lower 
stellar mass (for which our central values of ri are lower but 
poorly constrained ) sugg ests that our results and those of 
iRead fc TrenthamI ^00^ are consistent within the errors. 
Furthermore, their results and ours both suggest the trend of 
decreasing rj with increasing stellar mass above IO^/i-^Mq. 
They attribute this trend to feedback from AGNs which 
would te nd to lower the bar yonic fraction in high mass galax- 
ies re.g.. lSilk fc Reejll998D . 

4.5 Satellite fractions 

Here we present results for the satellite fractions as a func- 
tion of luminosity, stellar mass, and morphology. Table |S| 
and Fig. |7| show the best-fit a in stellar mass bins with 95 
per cent confidence intervals from the bootstrap, for early 
and late types separately. Table |7| shows the best-fit value 
of a in luminosity bins with 95 per cent confidence inter- 
vals, for early and late types, and for early-types split by 
local density, and Figure |H| shows these results as well. We 
note that Mcent and a have a cross-correlation coefficient of 
about —0.7 from the fits, because the choice of Mcent deter- 
mines the mass at which {N{M)) becomes proportional to 
Al, which has a significant effect on the satellite contribu- 
tion. 

We discuss both Figs. [7| and |H| together. For both the 
splits by luminosity and by stellar mass, the trend is for 
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Table 7. The satellite fraction determined for luminosity subsamples separately for early- and late-type galaxies, and for early-type 
galaxies split into low-density and high-density samples. All confidence intervals are 95 per cent. 
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Table 6. The satellite fraction determined for stellar mass sub- 
samples separately for early- and late- type galaxies. All confi- 
dence intervals arc 95 per cent. 
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Figure 7. Satellite fractions as a function of stellar masses for 
samples split into early (red hexagons, solid lines) vs. late (blue 
triangles, dotted lines) types. Errorbars shown are 95 per cent 
confidence intervals. 



a for early types to decrease slightly with L or M stellar, 
with central values around 0.5 at Msteiiar ~ 0.8 x IO^^Mq 
or L/L* ~ 0.1, decreasing to a central value of 0.1 (upper 
limit ~ 0.2) for M^tMar > 20 x IO^^A/q, or 0.2 (upper limit 
~ 0.35) for L/Lt > 2. Wo remind the reader that these 
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Figure 8. Satellite fractions as a function of luminosity for sam- 
ples split into early (red hexagons, long-dashed lines) vs. late 
(blue triangles, dotted lines) types, and early types split into high- 
density (magenta crosses, dashed lines) and low-density (black 
squares, solid lines) samples. Errorbars shown are 95 per cent 
confidence intervals. 



numbers do not limit the number of early-type lenses that 
are in groups and clusters, because a galaxy that is the BCG 
of a group/cluster (of which we expect many in L5 and L6) 
will not be included as a satellite. Hence, it is possible that 
the number of early-type galaxies in groups or clusters is 
actually constant with L or Msteiiar when the number that 
are BCGs is accounted for; unfortunately, we cannot test 
this hypothesis using our halo model formalism. 

For early types, we also have done a simultaneous split 
by luminosity and density. As shown, for L/L* <~ 2, this 
split seems to do an excellent job of isolating those early type 
galaxies that are in groups and clusters. For example, for L3 
(L/L, ~ 0.5), we see that the central value of a for the 
full sample of early type galaxies is 0.44, for the low-density 
elliptical sample is consistent with zero, and for the high- 
density elliptical sample is 0.96. Since the low- and high- 
density samples are determined by splitting at the median 
environment estimate, these results for satellite fractions are 
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fully consistent with each other. For all luminosity bins, sim- 
ilar consistency relationships are satisfied within the noise. 
We expect that for higher luminosities, many of the high- 
density sample are BCGs of groups and clusters, which may 
account for the discrepancy noted in central halo masses in 
the previous section, and for the decline of satellite fraction 
with luminosity for this sample. 

For late types, there is no suggestion of a statistically 
significant change in the group/cluster fraction with lumi- 
nosity or stellar mass, with central values ~ 0.10-0.15 and 
upper bounds typically 0.25-0.3. The lower satellite fraction 
for late types than for early types is consistent with works 
cited previously. 

The overall trend of satellite fraction decreasing with 
mass is also observed wit h simulations and semi-analytic 
galaxy formation models of IZheng et alj (|200^. 

We c ompare these r esults against derived satellite frac- 
tions from lZehavi et all (120051) for L2, L3, and L4. They find 
that 10-30 per cent of blue galaxies are satellites indepen- 
dent of luminosity, consistent with our results. Of the red 
galaxies, they find satellite fractions of 0.72, 0.54, 0.35 for 
L2, L3, and L4 respectively, so the trend of a decreasing 
with luminosity is found in both their autocorrelation anal- 
ysis and our lensing analysis. However, our values of a for 
these bins are 0.521^;^^, QAAt^fo, and 0.27l{];g5, about 30 
per cent lower than their results, though for L2 and L3 their 
values of a lie within our 95 per cent confidence intervals; 
for L4, they lie slightly outside our 95 per cent intervals, but 
when the errors on their estimates are included, it is not clear 
that the discrepancy is significant. They also find that the 
average host halo mass for L2 galaxies (i.e., the mass of the 
full cluster if they are in a cluster, or the halo mass of field 
galaxies, wh ich are in the minority) is2.5xlO"/i-^M0. We 
have found (iMandelbaum et al.ll2005bl) that the lensing sig- 
nal on group and cluster scales is actually dominated by less 
massive halos, in the mass range lO^'^-lO^^; i f the average 
host h alo mass was as high as the value given in lZehavi et alJ 
(|200^, then considering the large satellite fraction in this 
bin, the lensing signal would be too large to be consistent 
with observations on 500 /i^^kpc- 1 /i~^Mpc scales (increas- 
ing by more than a factor of two from its current value). 

4.6 Robustness of best-fit parameters 

In order to determine the robustness of these results, we 
consider which modeling assumptions may affect them. The 
first effect we consider is that of {N{M)) which has been 
modeled here as a power-law {N{M)) oc M for M > SMcent 
and oc for M < 3Mcent, so SMcent marks the point 
below which it falls off rapidly. This choice was justified in 
321 but we must consider the effects of this choice on our 
results for the best-fit satellite fractions. In particular, in 
[Mandelbaum ct al. (2005b), we found that the power- law 
exponent e of {N(M)) for M > SMcent was almost com- 
pletely degenerate with a. The two parameters essentially 
arranged themselves to preserve the amount of signal com- 
ing from halos in the mass range 10^"^-10^* /i~^Mq (higher 
mass halos are not important because they are killed off 
by the exponential in the mass function dn/dM, and lower 
mass halos do not contribute a significant amount of signal) . 
Thus, our satellite fractions have the potential to differ sig- 
nificantly from the real one if e is incorrect. 



Ivan den Bosch et al.l (l2003l) fi nd using d ata from 2 De- 
gree Field Galaxy Redshift Survey dColless et al. 2001, 2dF- 
GRS, ) that {N{AI)) has a different power-law dependence 
for spirals (shallower) and ellipticals (steeper), with the net 
result being that the bright luminosity bins, which are dom- 
inated by ellipticals, show a single power-law, but fainter 
luminosity bins are best described by a combination of the 
two power laws, with the spiral one dominating at low mass 
and the elliptical one dominating at high mass. This finding 
is reasonable in light of the fact that spirals are known to 
dominate in the field and ellipticals in clusters. However, it 
does mean that our assumption of e = 1 for both early- and 
late-types may complicate our analysis, since it may have 
caused an overestimate of a for early-types and underesti- 
mate for late-types, which is exactly in the direction of the 
difference we measured (i.e., lower satellite fractions for late 

types). 

As noted in Ivan den Bosch et al.l i2003i) , which seems 
to suggest e ~ 0.8 for ellipticals and ~ 0.6 for spirals using 
the 2dFGRS bj band data (in agreement with anot her bj 
band study using APM data. lScoccimarro et al .120011 . which 
found e ~ 0.8 overall), we do not expect e to be the same 
when determined using data selected by absolute luminosity 
in different bands, so in order to estimate its value for early- 
an d late-types and its p ossible luminosity evolution, we turn 
to IZehavi et al.l (1200 5^1. which uses SDSS samples selected 
in the r-band and a halo model analysis of ^gg{rp) in order 
to determine this parameter.^ The right-hand side of figure 
18 in that paper shows the best-fit value of e for samples 
selected by luminosity thresholds rather than 1-magnitude 
wide bins; as shown, the value of e is ~ 0.9 for Mr < — 18 
samples, rising slowly to e ~ 1.3 for Mr < —20.5 samples, 
then rising sharply to e ~ 2 above that. This trend may 
refiect the relatively higher fraction of early types in the 
brighter luminosity-threshold samples. The lower right panel 
of their figure 22 shows the best-fit {N{M)) for early- and 
late-types, from which we deduce values of 6 ■- ^ 1.65 and 
~ 1.10, respectively. We note that IZheng et alJ (|2o3), on 
the other hand, find e ~ 1 for all samples with the same 
data, reflecting a difference in th e modeling. The diff erence 
may be attributed to the fact that lZehavi et all (l2005l) lack a 
lower mass cut-off for their {N{M)}, whereas this work and 
IZheng et al.l (|20o3) do i nclude a lower mass cutoff. Because 
of the lack of cut-off in IZehavi et al.l (l2005l) . their value of 
e must necessarily be quite high for the brightest samples 
in order to l imit the contributio n from lower mass halos. 
We note that lConrov et aP (l2005h have explicitly tested this 
hypothesis by fitting for e with and without a lower mass 
cutoff, and found that the lack of lower mass cutoff can 
increase the best-fit e by as much as 50 per cent. 

We consider the effect of changing e on the best-fit val- 
ues of a for early and late types sepa rately by calculat- 
ing an analytic correction assuming the Zehavi et alJ (l2005h 
values of e for those samples independent of luminosity or 
stellar mass; this correction requires that we compute the 
non-central lensing signal for satellites residing in 10^"^-10^* 
h~^MQ halos as the product of a times an integral involv- 



1 Note that both lZehavi et alj i2005l) and lZheng et al.l J2004) call 
this parameter o (i.e., the symbol we use for the satellite fraction), 
rather than e. 
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ing {N{M)), and requiring that it be preserved when we 
change e, thus telling us the new value of a as well. We 
note that the correction factor is a function of central halo 
mass because, for high values of central halo mass, the cutoff 
SMcent will be within our mass range of interest, so changing 
e will have less of an effect. For both stellar mass and lumi- 
nosity bins, the results of this calculation indicate that for 
late types, ch anging e from our assumed value of 1.0 to the 
IZehavi et alj (^2005) value of 1.1 requires that we decrease 
a by 4 per cent of its best-fit value, which is significantly 
less than the statistical error on this value for any stellar 
mass bin. For stellar mass samples, for early types, chang- 
ing e from 1.0 to 1.65 requires that we decrease a by 21 
per cent of its original value for the 5 lowest stellar mass 
bins, 12 per cent for sm6, and no decrease for the largest 
stellar mass bin (due to its high central halo mass). For lu- 
minosity samples, for early types, it is to decrease a by 21 
per cent for the 5 lowest bins ( — 17 > Air > —21.5), by 
12 per cent for —21.5 > Mr > —22, and no decrease for 
—22 > Air > —22.5. Thus, the value of a for late types is 
nearly unchanged, and for early types, the limiting value at 
lower stellar mass becomes 0.4 rather than 0.5. We note that 
for the split by environment, since we only use early types, 
the same correction factors apply as for the full early type 
sample, so while the actual values of a must be lowered, 
this correction does not affect the consistency relationship 
between the satellite fractions in the low- and high-density 
early-type samples relative to the full early- t ype samples. 
Finall y, we remind the reader that since the IZehavi et alJ 
J2QQ^ values of e are higher than our assumed value in part 
due to a significant difference in modeling, and such high 
values of e are not consistent with N-body simulation re- 
sults when using a model similar to ours, these corrections 
are actually quite conservative. 

Next, we consider the fact that the expected signal de- 
pends on the radial distribution of satellites within groups 
and clusters. As mentioned in ^ this distribution is not well 
known, and is assumed here to be an NFW profile with the 
same concentration parameter as the DM. Variations in Cg 
cause the noncentral signal to peak a t differ ent characteris- 
tic radii as shown in lCxizik fc Selia!^ l)2002h . thus affecting 
primarily the shape of the signal rather than its amplitude. 
Unfortunately, since in most cases the values of a are ~ 0.2- 
0.3, we are not highly sensitive to Cg and cannot place much 
of a constraint on it. 

Several observational results have suggested that within 
clusters, re d galaxies are more centrally concentrated than 
blue ones llOemlei[ll974l : iMelnick fc Sargentlll977l : iDressleJ 
ll98nl : lAdami et al. ''l998Y Due to a relatively low sensitivity 
to Cg, we do not fit for it, but see what happens to the fit 
and to the best-fit Adcent and a if we use Cg — 2cdm ~ 24 
for early-types, and Cg — 0.5cdm ~ 6 for late-types. We do 
this comparison only for L3, L4, L5f, and L5b, since these 
are the samples with the greatest statistical power. 

For early types, we find that because increasing Cg 
moves the non-central contribution of the signal to rela- 
tively small scales, this change actually affects the central 
halo mass Adcent more than it affects the satellite fraction, 
with the tendency being to decrease the halo mass slightly to 
compensate for the higher non-central signal on 100 — 200 
/i~^kpc scales. This change only affects the best-fit by 
~ 1, with the change being in different directions for the dif- 



ferent luminosity bins. The masses in L3, L4, L5f, and L5b 
decrease to 60, 59, 87, and 86 per cent of their values from 
the fits with Cg = Cdm, and the satellite fraction changes by 
-0.01, -0.02, -0.01, and -0.01 (absolute value, not per cent of 
original). Hence, the changes to a for early types are well 
within the statistical errorbars, and we conclude that best- 
fit values are relatively robust to uncertainties in Cg. The 
changes to the best-fit values of Adcent are at most la (for 
L3) and usually somewhat less than that, and thus we con- 
clude that this parameter is also not sensitive to changes in 
Cg within our statistical errorbars, though with more data 
this statement may no longer be true. 

For late types, the decrease from Cg = Cdm to Cg = 
0.5cd77i decreased the best-fit by at most 1, and had less 
of an effect on the best-fit Alcmt because it shifted the non- 
central contribution to higher radius. The masses in L3, L4, 
L5f, and L5b increased by 9, 5, 1, and 15 per cent of their 
original values, well within the la errors, and a increased 
by 0.02, 0.02, 0.00, and 0.00, also weU within the errors. We 
thus conclude that uncertainty in the distribution of satel- 
lites within groups and clusters is not a significant source of 
systematic uncertainty in our estimates of satellite fractions 
and central halo masses. 

Another source of uncertainty in these estimates is our 
modeling of the h-h term, which affects the signal on large 
scales, and thus can change the best-fit a and (through their 
degeneracy) Alcent- Fortunately, since the h-h term is small 
on r < 2 h~^Mpc scales, we find that neglecting it entirely 
changes the best-fit parameters by less than la. 

Finally, we remind the reader that due to the assump- 
tion of a universal IMF when deriving stellar mass estimates, 
these estimates may need to be rescaled by a constant fac- 
tor of up to 30 per cent (conservatively) . The main results 
that are affected by this rescaling are the limits that we have 
placed on conversion efficiencies. All other mass-related re- 
sults, such as trends in halo masses with stellar mass, and 
comparisons between different morphology or density sam- 
ples at constant stellar mass, are unaffected. 



5 CONCLUSIONS 

In this work, we have used halo model analysis of the galaxy- 
galaxy weak lensing signal in order to observe trends in the 
relationship between stellar masses and halo masses, and 
luminosity and halo mass, treating samples based on mor- 
phology separately. We have also studied ellipticals in low- 
and high-density regions separately. 

As a result, we have come to a number of conclusions 
related to average halo masses. First, the AI stellar /Al cent ra- 
tio is highest for AAgteiiar ~ 10^^ Mq, with a peak conversion 
efficiency of roughly 14 per cent for ellipticals (> 11 per cent 
at 95 per cent CL, statistical error only), although we can- 
not exclude even larger values at low luminosity or stellar 
mass where the lensing signal is weak. The corresponding 
number for spirals reaches 35 per cent at the maximum, but 
with a larger measurement error, which allows it to be as 
low as 18 per cent (95 per cent CL). These results imply 
a factor of two or more difference in conversion efficiency 
between typical spirals and ellipticals above stellar mass of 
10^^ Mq, whereas below this stellar mass, we find no incon- 
sistency between the conversion efficiencies of early and late 
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types, implying that stellar mass is a good tracer of halo 
mass in this regime. Note that since we are only including 
central galaxy stallar mass in the analysis, our conversion 
efficiencies should be viewed as lower limits, and also serve 
as lower limits on star formation efficiency. In practice, for a 
10^^ Mq halo, the central galaxy likely contains most of the 
stars and the contribution from satellites is negligible, while 
for a 10^^ Mq halo other galaxies outside the halo center 
contribute significantly to the total stellar content of a clus- 
ter. Thus while the upper limits to the conversion efficiency 
are uncertain both because of large observational errors and 
because of modeling uncertainties, a lower limit of about 
10 per cent (95 per cent CL) is robust for both spirals and 
ellipticals. Similarly, M/ L reaches a minimum for L Lt 
galaxies of 41 j^ij?-^©/^© for late types, or for L ~ L,/2 of 
51I37 for early types (95 per cent CL). Below these values 
of stellar mass or luminosity, both of those quantities trace 
halo mass well, meaning that halo masses for early and late 
types were consistent. At higher stellar mass or luminosity, 
the early types have a larger central halo mass, likely refiect- 
ing the fact that they tend to reside at the center of clusters, 
unlike late types. 

We also have a number of conclusions regarding satel- 
lite fractions. The satellite fraction for late types tends to 
be 10-15 per cent regardless of luminosity. For early types, 
the situation is more complex, with a satellite fraction of 
^ 40 — 50 per cent at low luminosity, decreasing to about 
20 per cent at high luminosity. These trends are consistent 
with those based on auto-correlation analysis iZehavi et alJ 
[200^, though the actual values of satellite fraction for early 
types are slightly different, possibly reflecting differences in 
modeling. Our fit results for L <'~ L, early types indicate 
that our density estimate is highly efficient at isolating a 
nearly pure satellite sample in this regime, which can then 
be used in a future work that will study tidal stripping and 
the radial distribution of satellites within groups and clus- 
ters. Current results suggest that tidal stripping is not com- 
pletely efficient in removing the dark matter from the satel- 
lites, since the satellite sample shows plenty of lensing signal 
at small (50-100 /i~^kpc) scales. 
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